Propulsion of a molecular machine by asymmetric distribution of reaction— products 
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A simple model for the reaction-driven propulsion of a small device is proposed as a model for 
(part of) a molecular machine in aqueous media. Motion of the device is driven by an asymmetric 
distribution of reaction products. The propulsive velocity of the device is calculated as well as the 
scale of the velocity fluctuations. The effects of hydrodynamic flow as well as a number of different 
scenarios for the kinetics of the reaction are addressed. 
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Molecular motors are machines that convert chemi- 
cal energy to mechanical work [J Q. Examples are 
the cytoplasmic motors that move along biological (pro- 
tein) tracks in the cell by converting the energy re- 
leased upon ATP hydrolysis into mechanical work [3, 0] • 
These complex machines act as the inspiration for the 
design of macromolecular devices [j| with the ability to 
sort, sense and transport material in chip-sized laborato- 
ries [6(1 . Consequently a major area of current chemical 
research is the construction of much simpler molecular 
machines for nanotechnological applications Q- 

In this spirit, we study in this Letter a simple model 
of a self-propelling device driven by chemical reactions 
on its surface. It is simple enough that not only is the 
construction using chemical techniques feasible, but it 
should also be possible to change parameters in order to 
optimize particular features or functions. 

We consider a spherical particle (colloid or vesicle) of 
radius R that has a single enzymatic site located on its 
surface at a fixed position, as sketched in Fig. [TJ In 
the presence of a reactive substrate in a non-equilibrium 
state, the enzyme promotes the reaction rate in its vicin- 
ity and produces a dynamic and asymmetric distribution 
of product particles of (hydrodynamic) diameter a«i? 
which exert osmotic or interfacial forces (depending the 
boundary properties) and hence propel the sphere in a 
fixed direction. We consider both uniform and peri- 
odic reactions and calculate the propulsive velocity of 
the sphere. We find that the velocity of propulsion is 
set by the size of the product particles, the properties on 
the boundary, and the reaction rate. Upon considering 
the variations in the rate particle release (and taking ac- 
count of density fluctuations), we find that the ratio of 
the mean-square velocity fluctuations to the mean veloc- 
ity depends on the ratio of the time for a product particle 
to diffuse a distance R to the typical time between the 
production of successive product particles. 

The reaction site on the sphere, located at r s = — iR 
(see Fig. [IJ, is an enzyme catalyzing the break up of 
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FIG. 1: The spherical particle with a reaction site. 



an available substrate into two product particles. If one 
of the product particles is similar in size to the original 
substrate (see below), we can assume that the reaction 
site effectively acts as a source of excess particles that 
are being released at the reaction site at a rate dN ?} , 
and the diffusion equation for the density of the excess 
product particles can be written as 



d tP (r,t)-DV 2 p(r,t) 



dN p (t) 
alt 



£ 3 (r-r s ), (1) 



where D is the diffusion constant of these particles. The 
density profile can be obtained from Eq. (JTJ) subject to 
the boundary condition of vanishing normal current on 
the surface of the sphere, f • Vp| sp horc = 0. 

The resulting distribution of product particles around 
the sphere is asymmetric, with a non-zero first moment 
of pi(i) = J sin 6 d8 cos9 p(r = R,6,t), leading to 
phoretic motion of the sphere. It is well known that 
colloidal particles in externally imposed solute gradients 
will be set in motion by a variety of phoretic mecha- 
nisms ||. Here in contrast, the gradient (of products) 
is self- generated by the device itself. We thus ob- 
tain a general expression for the velocity of propulsion in 
the z-direction v = — ^ p%, where n is the viscos- 
ity of the solvent and the length-scale A depends on the 
particular phoretic mechanismj^, |T(|. The other com- 
ponents of the velocity vanish due to symmetry. Two 
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mechanisms valid for non-ionic product particles are the 
"diffusiophoresis" of a totally impermeable sphere 



because the depletion of the product particles near its 
surface causes a lateral slip velocity that results in net 
motion of the sphere, and the "osmiophoresis" of a spher- 
ical shell which is permeable to solvent but impermeable 
to product particles which develops a non-zero velocity 
due to osmotic forces that cause radial flows of solvent 



across the membrane[12l|. 

Experimental estimates for the diffusiophorctic, Ao 
and the osmiophoretic, Ao lengthscales are available from 
the literature. Measurements on latex particles in gradi- 
ents of Dextran by Staffeld and Quinn [13] obtained val- 
ues of Xd — 38 nm. Studies of micron sized lipid vesicles 
of radius R in gradients of sucrose by Nardi et al Q have 
obtained an osmiophoretic lengthscale Ao ~ 0.3R. Theo- 
retical treatments of phoretic propulsion implement local 
momentum conservation by solving the Stokes' equation 
for the solvent (Re<§; 1) taking account of the product 
particles and force balance on the sphere 0. For dif- 
fusiophoresis, Anderson and Prieve [l5| have calculated 
X% = J °° dll [1 - e - w ^l k * T ] , where W{1) is the inter- 
action energy between the solute particles and the rigid 
wall of the sphere at a normal separation I This is 
consistent with the experiments above [l3T | for an inter- 
action range given by the hydrodynamic radius of Dex- 
tran. For osmiophoresis, Anderson [l§] has obtained, 



R 2 



r/Lp/R 



where L p is the filtration coef- 



For a R = 2^m lipid vesicle, 
10~ 7 (m/s)/atm [16J, one 
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ficient of the membrane 
using i] = 10~ 3 Pa.s and L 
obtains ~ 0.01 significantly lower than the experi- 
ments above [3]. This discrepancy between theory and 
experiment remains an open problem [lbT |. 

Solving Eq. ([T|) subject to the appropriate boundary 
condition mentioned above, we find p±, which together 
with the above expression for the velocity of propulsion 
yields 



v(t) = I 



d'jj 



7M 



jR 2 



2tt 



^R 2 
D 



^oR 2 
D 



(2) 



where £ — and f(u) is the Fourier transform of 

dN j} ■ Using Einstein's formula for the diffusion coeffi- 
cient D — kftT I (Sirrja) , we find 



1 = 



3 A 2 
4i? 2 



(3) 



The time dependence of the velocity will depend on the 
particle release rate. Two important cases of uniform and 
periodic particle release are considered. 

Uniform Particle Release. Consider the case where 
a particle source with a uniform rate is switched on at 
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FIG. 2: Transient response of the velocity of propulsion after 
the reaction has been switched on. 



t = 0; i.e. — where T{ is the average in- 

verse reaction rate or "firing time" of the product par- 
ticles, and d{t) is the Heaviside step function. Using 
Eq. ([2]), we find the average velocity of the sphere as 
v(t) = v Q (Jp) , where the stationary velocity is 
given by 



"o = — , 



(4) 
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on [13 . The 



cos 2s -sin 2s + 4^ iF 2 (1; |, |; -s 2 ) , 
s 2 ) is a generalized hypergeometric 



and Q(s) 
where 

function [17] . The time response of the force is plotted 
in Fig. The function starts at t = with an infinite 
slope and asymptotes to its final value around t = R 2 /D. 

We can estimate the steady state propulsive velocity 
for a spherical device of radius R = 2/xm with the en- 
zyme on its' surface catalyzing a fast reaction such as that 
of acetylcholinesterase, which hydrolyzes acetylcholine in 
the synaptic cleft into components acetate and choline 
at the rate of 1/rf ~ 25000 s _1 (HI- The corresponding 
product particle will be acetate with a hydrated diam- 
eter of a = 0.8 nm. Using the estimates above 13|, 15 1 
for Afl we find a steady-state propulsive velocity from 
diffusiophoresis of a hard sphere as vu ~ 1 nm/s. This 
rather small value could be improved by "engineering" 
the strength and range of the interaction of products 
with the surface. The high experimental value of Ao 
measured [3] for closed vesicles imply a much higher ve- 
locity of propulsion vo ~ 1 um/s for a semipermeable 
shell. Theoretical estimates [12| however predict a much 
smaller velocity vo ~ 1 nm/s. 

Periodic Particle Release. If the reactions oc- 
at well defined time intervals of Tf, then 

velocity reads v(t) = 

" 1 ■ ' (*) 



cur 




1 + (1 — i)y/Timx — iirmx 
be simplified in two limiting cases 



where A r 



The above expression can 



adt7 



<C 1 , the ex- 



pression in the brackets can be approximated as 1 and the 
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FIG. 3: Time evolution of the propulsive velocity for the case 
of periodic particle release, corresponding to R/\JDrt = 1. 

series can be summed up to v(t) ~ v J2n=-oo <KV r f — n )- 
This corresponds to the case where the diffusion time 
around the sphere is much less than the firing time. In 
this case we expect the initial inhomogeneous profile of 
the product particles to lead to instantaneous impulses 
that are immediately screened due to fast homogeniza- 
tion. In the opposite limit of 3> 1, however, the 

rapid release of the product particles maintains a stabi- 
lized inhomogeneous profile and the velocity is approxi- 
mately constant: v(t) ~ vq. The time evolution of the 
velocity is plotted in Fig. [3] for an intermediate case of 
—M= — 1. We note the that the relative magnitude of 
the velocity fluctuations and the average velocity is set by 
the ratio between the firing time and the diffusion time. 

Velocity Fluctuations. The velocity calculated in Eq. 
(|2|) should be considered as the average of a fluctuat- 
ing quantity. There are two sources of fluctuations that 
need be taken into consideration for a quantitative assess- 
ment of the velocity fluctuations, namely, product parti- 
cle density fluctuations and randomness in the reaction 
that leads to the particle release. 

To take account of the density fluctuations, we can go 
back to the Langevin dynamics of the particles described 
as dtYi(t) = u.i(t) for the velocity of the i-th. particle, 
where u^(t) is a random noise with a distribution P[u] = 
cxp [— Yli I dt Ui(i) 2 ] . We can construct a stochastic 
density as p(r,i) = X^ 3 ( r — r -»(*))j an< ^ snow that it 
satisfies Eq. ([!} with an additional noise term Q(r,t) 
added to its right hand side, whose moments can be 
calculated using the above distribution as (Q(r,t)) = 
and (Q(r,i)Q(r',t')> = 2D(-V 2 )S 3 (r~ r')S(t- t')p(r, t), 
where p(r,i) = (p(r,t)). To incorporate the random- 
ness of the reaction, we write the particle release rate as 

dN ^ l) = J2 n <K* _ Hj<n r i)' wnere we have defined a 
time interval r„ between the release of the n + 1-th and 
the n-th particles. We further assume that the reaction 
leading to the product release is a Poisson process, in 
which case the probability distribution of the time in- 
tervals between two consecutive particle release events is 
given by P{t) = — er T l Ti . 

Using the formulation outlined above (details of which 
will be presented elsewhere [loj). we can calculate the 



fluctuations of the velocity about its mean value. We 
find the long-time behavior of the velocity correlations 
as ([v(t) — v ][v(t') — v ]) ~ ^ S(t — t r ), where the over- 
line denotes averaging with respect to the time-interval 
distribution. This result implies diffusive behavior, and 
we obtain an effective diffusion coefficient T> zz ~ ^ 2 /t[, 
for the sphere in the propulsion direction. Similar re- 
sults are obtained for the effective diffusion coefficients 
V xx and Vyy in the lateral directions. We can also es- 
timate the relative importance of the velocity fluctua- 
tions as compared to the average velocity by calculating 
(Aw) 2 = ([«(<) — vq} 2 ). In the long-time limit, we find 
^ Av } ~ Dt(/R 2 , for R ^s> a. The velocity fluctuations 
are thus controlled by the ratio between the firing time 
and the time it takes for the product particles to diffuse 
across the sphere and homogenize their profile. For ro- 
bust propulsion, the firing time must be considerably less 
than the diffusion time. 

Effect of Hydrodynamic Flow. To make the calcula- 
tions self-consistent we should also take account of the 
fact that the diffusion of the product particles will be af- 
fected by the hydrodynamic flow around the sphere. The 
velocity profile of the solvent around the sphere can be 
shown to have a 1/r 3 dependence, which is character- 
istic of phoretic flows (l3 |. A corresponding convective 
term should then be added to Eq. {T]), which suggests 
a systematic solution of the equations as a perturbative 
series in the Peclet number Pe = vR/ D. Previous work 
has shown that the first correction term is of the order 
of Pe 2 for osmiophoresis [T^, and Pe (Xd/R) for diffu- 
siophoresis [15}. These corrections are negligible for the 
cases considered here (see below). 

Reaction Kinetics. It is illuminating to examine more 
closely the conditions under which the above simplified 
picture holds, taking into account the reaction kinetics. 
We have a solution of substrate S at a relatively high 
concentration Cs that has a natural slow tendency to 
dissociate, S — ► S' + V where S' and V represent the 
product particles that are assumed to exist in vanish- 
ingly small concentrations in the bulk. This means that 
the entire bulk of the system is in a non-equilibrium con- 
dition that will not relax to equilibrium in laboratory 
time scale. The presence of an enzyme E catalyzes this 
reaction, so that it takes place considerably faster in the 
vicinity of the enzyme, which leads to an increased pro- 
duction of iS' and V with corresponding concentrations, 
Cs 1 , and C-p respectively. There will also be a depletion 
of S in the neighborhood of the enzyme. The products 
will slow down the reaction by moving it towards equi- 
librium using the backward path S' + V — ► S that is 
inevitably present, but one can show that it is only a 
small perturbation to the reaction condition, because all 
the concentrations have to match smoothly with those of 
the bulk that are maintained out of equilibrium. The net 
propulsion of the sphere depends on all three concentra- 



4 



tions. If the diffusion constants (i.e. the hydrodynamic 
radii) and other properties of the S and the S' particles 
are not that different, the sum of the two concentrations 
Cs + Cs> remains constant and their contributions to the 
phoretic propulsion cancel . This leads to the simpli- 
fied picture presented above of a diffusive process with a 
localized source where p — C-p in Eq. {T]). 

To gain further insight on what controls the effective 
particle release rate, we consider the following general 
multi-stage reaction pathway 



S + £ 



S£ 



k 2 



V 2 £ + Vi^Vi+V 2 +£, (5) 



where the two products Vi (1st) and V 2 (2nd) repre- 
sent S' and V , the order depending on which of them 
{S 1 or V) is released first. The backward reaction paths 
in Eq. ([5]) have been eliminated for simplicity, and 
they can be shown to have negligible effect provided 
the noncquilibrium working condition described above is 
maintained [1 01 ] - Under steady-state reaction conditions, 
the diffusion-reaction equations take on the form of Eq. 
(fT]), with the same rate (as a sink for the substrate and 
source for the products) appearing on the right hand side 
of these equations. The mean rate constant is given by 
the Michaelis-Menten rule [2(1 



f 

Tf 



Cs{£) 



k 2 + k 3 J Ku + C s (Sy 



(6) 



where K M = ^ 



k2 + ks 



is the Michaelis constant and 



Cg{£) is the substrate concentration at the position of 
the enzyme. 

The concentration of the substrate at the position of 
the enzyme, which we expect to be depleted in compari- 
son to the bulk concentration Co, can be found by solving 
the corresponding reaction-diffusion equation with the 
appropriate boundary condition. We find 



Cs(£) = = 
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where as is the hydrodynamic radius of the substrate 
and ag is a typical size of the enzyme. 

Discussion. The molecular locomotive device will be 
useful if it can perform directed motion over a distance 
that is large compared to its own size. A very im- 
portant limiting factor for directional motion is rota- 
tional diffusion, whose time scale for a spherical colloid 
is tr = SirnR 3 / (kftT), which implies that we can achieve 
more directed motion by increasing the size of the sphere. 
The rotational diffusion time turns out to be of the order 
of 50 s, for R = 2 fim, which means that given an optimal 
propulsive velocity of vq ~ 0.5 /im/s, a locomotive of that 



size can travel of the order of ten times its own size be- 
fore it loses sense of its original orientation. The effective 
diffusion constants derived from the velocity fluctuations 
above contributes negligibly to the orientational mem- 
ory loss, giving a rotational diffusion time of the order of 
t v f ~ (R/£) 2 Tf ~ 1000 s. We also note that the typical 
time that the system spends in a transient regime after 
a change is given by tq = R 2 /D (see Fig. [2]), which is of 
the order of 5 ms. The effect of the hydrodynamic flow 
can be seen to be unimportant, because Pe ~ 10~ 2 . 

In conclusion, we have proposed a model design for 
a molecular machine that can propel a vessel in aque- 
ous media with a mechanism that involves asymmet- 
ric release of reaction products. The motor may be 
thought of as a diffusive equivalent of the jet engine: it 
releases asymmetrically the reaction products in a vis- 
cous medium, lets them diffuse and takes advantage of 
their thermodynamic forces, instead of gaining inertial 
thrust by ejecting the exhaust. 
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